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Abstract 



We propose a new method to investigate collective behavior in a network 

of globally coupled chaotic elements generated by a tent map. In the limit 

of large system size, the dynamics is described with the nonlinear Frobenius- 

Perron equation. This equation can be transformed into a simple form by 

making use of the piecewise linear nature of the individual map. Our method 

is applied successfully to the analyses of stability of collective stationary states 

and their bifurcations. 
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I. INTRODUCTION 



Globally coupled dynamical systems form an important class of models in nonlinear 
dynamics. The all-to-all nature of the coupling may be regarded as an idealization of long- 
range coupling or an approximation to a short-range coupling high-dimensional lattice. Such 
systems are useful in modeling diverse phenomena such as Josephson junction arrays mul- 
timode laser [Q, charge- density wave evolutionary dynamics 0], biological information 
processing and neurodynamics A simplifying assumption is often made that the cou- 

pling strength is uniform over all coupled pairs. This makes the theoretical treatment much 
easier compared with systems with short-range interaction [0]. 



In the present paper, we study globally coupled maps (GCM) introduced by Kaneko [jT2 
They form a dynamical system of N local mappings which are under a common internal 
field. The system is thus a mean-field version of coupled map lattices (CML) The 
explicit form of GCM we will work with is given by 

N 



Xn+i{t) = (1 - e)/(x„«) + - ^ (1) 



e 

r 

i'=l 

Here n represents a discrete time step, i the index of the elements {i = 1, 2, 3, A^), and e 
the coupling constant. It is assumed that the map f{x) gives rise to chaotic dynamics. We 
specifically consider a tent map 

f{x) = l-a\x\ (2) 

with parameter a. The piesewise linearity of the map simplifies the analysis considerably as 
we see later. The range of x of each map is given by [—1, 1]. The coupling is only through 
the mean field 

1 ^ 

hn = ^J:fM^)). (3) 

i=i 

There are a great variety of phenomena generated by such globally coupled maps, and 
they seem to be due to two conflicting natures involved in the dynamics |]T2[. On one 
side, the presence of a common driving force coming from the mean field favors mutual 



synchronization. On the other side, unstable growth of the difference in x„ between any 
pair of elements tends to destroy such synchrony. If the coupling constant e is sufficiently 
large, mutual synchronization will be complete, all elements behaving identically. Then 
the collective motion is described simply by the single map f{x). On the contrary, if e is 
sufficiently small, the system will be "turbulent" where the elements are scattered and behave 
chaotically in time. In the present model there also appear band structures for intermediate 
coupling strengths. The state of perfect coherence loses stability at a(l — e) = 2, while the 
transition between turbulent and two-band states occurs at a(l — e) = V2. A phase diagram 



was obtained by Kaneko [T^ 



In this paper, we shall confine our analysis to turbulent and two-band states. The limit 
of large is of our particular concern. Some unexpected features of collective quantities 
such as the mean field have been found recently. For instance, the mean square deviation 
of such quantities dose not tend to zero as — oo but saturates to a finite value. This 
suggests the existence of nontrivial collective behavior p!3|-pO|]. 



II. NONLINEAR FROBENIUS-PERRON EQUATION 

For a very large system size, direct numerical simulation of (|l]) requires enormous com- 
puter capacity. Moreover, the occurrence of long transients may complicate the analysis of 
stationary states and other asymptotic states. For these reasons, an alternative approach 
seems desirable. 

Statistical states of an ensemble of maps governed by (|1|) can be characterized by its 
density distribution 

N 

N 

where 6{x) is Dirac's delta function. The evolution of pn{x) obeys the Frobenius- Perron 
equation ||T^p!5 



Pn(2;) = ^^5(a; - x„(i)), (4) 

i=l 



pn+i{x)= 5{x-Fn{y))pn{y)dy, (5) 



where the effective map F„(x) depends on the mean field hn through 

Fn{x) = [1 - e){l - a\x\) + eK. (6) 
The mean field hn is calculated from the density as 

hn= f{x)pn{x)dx. (7) 



Equation is called the nonlinear Frobenius- Perron equation (NFPE) In the limit 

of large system size the density (|) becomes a "sufficient continuous" function. NFPE is 
completely analogous to the nonlinear Fokker-Planck equation. In the following, we assume 
that the initial condition is given by a uniform distribution Pq{x) = |, that is, the elements 
of the original GCM are distributed in a completely random fashion. 

Owing to the linearity of the two branches of the tent maps, is simplified as 

Pn+l{X - 7- 7 , [H) 

a(l — e) 

where ?/+ and y_ are the two preimages of the same x, i.e. x = Fn{y±). 

In the following, we propose an expansion form of NFPE in order to facilitate the analysis 
of the model. Figure |I| shows snapshots of Pn{,x) for two parameter conditions. The density 
distributions for a turbulent state and a two-band state are displayed in Figs.[l|a and ||b, 
respectively. We assume that the density distribution is a piecewise constant function, 
which is due to the piecewise linearity of the tent map. Then, one may simplify (H) by using 



a technique adopted to analyze a one-dimensional tent map ||2T]|. The density distribution 
is expressed as 



where 6{x) is a step function 



Pn{x)=Y.KO{<-x). (9) 



1 (x > 0) 

e{x) = { - (10) 

(a; < 0) 
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In the above equation, 6^ and c{ represent the height and the location of the jumps in Pnix), 
respectively. Since the initial condition for the density distribution is assumed uniform, lf{ 
and Ci are given by 

6? = ^ 



a(l - e) 
ae 



others 



1 



a(l — e 
ae 



(11) 



ail-e) 







For n = 1, the number of jumps in the density Pn{x), which is denoted by equals two, 
and increases with n as u = n + 1. With the use of (|), the evolutional equations for and 
are given by 

2 

=0 



J+1 



a{l - e) p 
Fn{0) 

1 

~a{l-e' 



(20(c^j - 1) 



(12) 



This set of equations gives an alternative representation of NFPE (^. Note the identity 



i=o 



(13) 



which results from the requirement p„ 
h„ is calculated as 



0. By using (HI) along with (01), the mean field 



(14) 



j=0 



Equations (^, (|T2]) and (|T^ constitute a closed set of exact nonlinear equations equivalent 
to NFPE, and (p^l]) specifies the initial condition. 

The dimension of the dynamical system (|T2p is potentially infinite {j = 0,1,..., cxd) 
because the number of jumps in p„(x) increases indefinitely with n. We now approximate 



this by a finite number. We expect that the jumps of sufficiently large j hardly contribute to 
the collective behavior, so that we may ignore such minor jumps. Indeed, the third equation 
in (p!2| ) means that the height of the jumps is given by 

\K\ = N- lCl- (15) 

Thus, 6{ approaches exponentially as j oo, provided a(l — e) > 1. We now focus on 
this incoherent regime . Note that in the opposite case, i.e. a(l — e) < 1, p„(x) becomes 
a delta function since all elements are in perfect synchrony, so that the breakdown of the 
above approximation is apparent. Our approximate method seems more convenient than 
those employed previously [0,0 • In what follows, we restrict the maximum number of the 
jumps to 50 The mean field behavior calculated with this approximation turns out 

consistent with our direct numerical simulation. 



III. STABILITY OF STATIONARY STATES 

We investigate stationary solutions i.e. fixed points of NFPE which are however not the 
fixed points of GCM The fixed points of NFPE are determined by pn+i{x) = pn{x) = 
p^{x). From the expansion (|12D in the last section, this condition is given by 
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6? 



K*' = —r, -Jf.m'i)-l), 

a(l — e) 

where the map evaluated at the fixed point is given as 

F,{x) = {I - e){l - a\x\) + eK, (17) 
h. = {l-^Y.V.<W.\). (18) 
The above equations constitute a set of self-consistent equations. For the time-being, we 



regard h^, in ([T7|) as an input field and distinguish it from in (0) which we regard 



as an output field hout- We can easily calculate bl and ci from ([I6D and the normalization 
condition for the density 



^6i(ci + l) = l. 

j=0 

With bi and ci thus obtained, the static field hout is given by 



(19) 

). Consequently, hout is 



obtained as a function of i.e. hout{hin]a,e)^ which is called a "static mapping" [jT9 
The fixed points of this mapping give the stationary solutions of our system. The mapping 
actually has some fixed points in most of the parameter region under consideration. 

Numerical stability analysis of the stationary states is now presented. We first linearize 
(|l^) about a stationary state (6^, c^) corresponding to the density p*(a:), and consider a 
small density deviation from it. Noting that Fn{x) in (|12|) depends on 6{ and c{, we have a 

set of linear maps given by 

2 



ad - el 



7°-M = -f E{/?;i'civ:i+27;iv,v:i} 

^ i'=o 



i'=o 



Pn+l 



i+1 
Tn+l 



1 



(20) 



a(l - e 
ae 



■/?^(2^(ci)-i: 



. E { ci' Ici'l + 27^ |ci'|} - a(l - 6) 7^ (2^(ci) - 1) 
i'=o 



where /3i = - bl and 7^ 
coefficient matrix D as 



c^. The set of equations above can be expressed with a 



(21) 



The eigenvalues of D determines the stability of the stationary solution. It is obvious from 
fEOl) that D has the form 







= D 








y7nH 






[in) 



D 



V 



(22) 



* D2 

Thus, the eigenvalues of D are given by those of Di and D2, and they are generally complex 
numbers. In Fig.|^, we show the eigenvalues of D for a = 1.5 and e = 0.05 for which a 
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stationary state exists stably. They include real numbers and complex-conjugate pairs, and 
their distribution is roughly given by two concentric circles. The inner and outer circles 
correspond to Di and D2, respectively. The normalization condition for p„(x) requires that 
one of the eigenvalues is always unity [0]. If all of the remaining eigenvalues lie within 
the unit circle, the stationary solution remains stable. It loses stability in the following two 
ways. 

The first case is due to band splitting. This occurs at a(l — e) = |2^. In the unstable 
region, period-two motion appears, and then the density pn{x) has a two-band structure as 
shown in Fig.|I|b. An infinite number of period-two solutions exist and a particular state 
is selected out by imposing a specific initial condition. Therefore, it is expected that their 
stability is marginal. 

Figure ^ shows the eigenvalues under the condition that stationary state is stable (a(l — 
e) < V^). The moduli of all eigenvalues except for one with vanishing argument are smaller 
than unity. Figure shows the eigenvalues after band splitting (a(l — e) > V2). The matrix 
Di has two eigenvalues with moduli 1, while their arguments are given by to and vr (the 
eigenvalues are 1 and —1, respectively). In Fig.||, we plotted the eigenvalue with argument 
vr as a function of e. As is seen from Fig.^a, \X-„\ approaches unity as e increases up to the 
band-splitting point eo under fixed a at 1.5. For this value of a, the band splitting is seen for 
e > eq = = 0.05719 ■ ■ ■, and lA^rl is kept at unity over this range of e. The results of 

the stability analysis agree with our observation from direct numerical simulation. In order 
to investigate the nature of this instability in some detail, we make a log-log plot near the 
band-splitting point. This is shown in Fig.pD. As is seen from Fig.^D, |A,r| approaches unity 
linearly with eo — e. Thus, the band-splitting point can be determined accurately from this 
kind of analysis. 

There is another type of instability. This is associated with a Hopf bifurcation, and 
occurs even in the region where the band splitting is not seen (a(l — e) < a/2). Figure ^ 
shows the distribution of the eigenvalues in the unstable regime. IAt^I is smaller than unity 
and the matrix D2 has an eigenvalue with modulus larger than unity. The argument of 
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this eigenvalue is about O.SGtt. In this region, the collective behavior of the original system 
is found quasiperiodic. The evolution of the corresponding mean field /i„ is displayed in 
Fig.^, and its power spectrum is given in Fig.^. The latter has some sharp peaks, and the 
frequency of the highest peak is about 0.43. This value is consistent with the results of our 
stability analysis. 

IV. SUMMARY 

The piecewise-linearity of the map allows us to make an expansion of NFPE in term of 
the height and location of the jumps in the density distribution. We analyzed the stability of 
stationary solutions by using this expansion. Linearization of the expansion is expressed with 
two matrices Di and The band-splitting instability is associated with the eigenvalues of 
Di. Stability analysis predicts the band-splitting point accurately. The band splitting results 
in marginally stable period-two states. There is also a region where quasi-periodic motions 
appear. The corresponding instability is essentially a Hopf bifurcation, and associated with 
the eigenvalues of D2. 

Our analysis presupposes infinite system size. In actual finite system, however, we expect 
some fluctuations around the solutions of NFPE. 

We investigated only globally coupled tent maps. It should be mentioned, however, that 
the same methods could also be useful for investigating other piecewise linear maps. 
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FIGURES 

FIG. 1. Snapshot of the density distribution Pn{x) after transients have subsided, (a) Turbu- 
lent state obtained from the expansion of NPFE (see Section 2); a = 1.6, e = 0.05. (b) Two-band 
state; a = 1.6, e = 0.15. 

FIG. 2. Eigenvalues of the 50 x 50 matrices Di and D2 for a = 1.5 and e = 0.05 where the 
stationary solution remains stable. 

FIG. 3. Eigenvalues of the matrices Di and D2. Their absolute values vs. arguments are 
shown. Since the eigenvalues appear as complex-conjugate pairs, only those with arguments lying 
between and (a) a = 1.5, e = 0.05. The moduli of all eigenvalues are smaller than 1 except one 
with vanishing argument, (b) a = 1.5, e = 0.06. There exists an eigenvalue of modulus 1 and 
argument tt. (c) a = 1.9, e = 0.12. There exists an eigenvalue with modulus is larger than 1. 

FIG. 4. (a) Eigenvalue lA^rl with argument tt vs. e where a = 1.5. (b) Log-log plot near the 
band-splitting point eo where a = 1.5. 

FIG. 5. Trajectory of the mean field /i„ obtained from the expansion of NPFE (see section 2). 
The parameter values are the same as in Fig.2c (i.e. a = 1.9, e = 0.12). This was obtained from 
1000 steps after transients have subsided. 

FIG. 6. Power spectrum of the time series of the mean field hn obtained from the expansion 
of NPFE (see section 2). This was obtained from an average over 10 runs each of 1024 iterations 
after transients have subsided. Parameter values are the same as in Fig.2c. 
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